Data fusion using consensus aggregation functions

ABSTRACT

A fusion system fuses M rankings generated by M judges by (i) computing values of an aggregation function for items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item, and (ii) constructing an aggregation ranking based on the aggregation function values. In an illustrative application, the judges are different Internet search engines and the rankings are sets of search engine results generated for a query input to the search engines, and a consensus search result corresponding to the query is defined by the aggregation ranking. In another illustrative application, the judges are different soft classifiers, and the rankings are probability vectors generated for an input object by the different soft classifiers, and the input object is classified based on a consensus probability vector defined by the aggregation ranking.

BACKGROUND

The following relates to the informational arts, computer arts, and related arts. Some illustrative applications of the following include meta-search engines, meta-classifiers, meta-prioritizers, and so forth.

In the informational arts, a common operation relates to scoring items based on a scoring criterion. For example, an Internet search engine receives a scoring criterion in the form of a search query, and outputs a list of relevant web sites coupled with relevance rankings or scores, where the relevance ranking or score assigned to each web site reflects the search engine's assessment of the relevance of that web site to the search query. In a simple embodiment, the relevance ranking or score may be a count of the number of occurrences of terms of the search query in the content of the web site. More complex relevance scoring or ranking may take into account other factors such as where the search terms occur (for example, search terms in the title of the web site may be weighted more strongly than search terms in body text), the popularity of the web site (measured, for example, based on a count of “hits” or downloads of the web site), or whether the web site is a “favored” web site (for example, a web site that has contracted with the search engine operator to obtain more favorable relevance rankings). In general, different search engines may use different bases for computing the relevance.

Other examples of scoring items based on a scoring criterion include: soft classification in which a classifier assigns an object a probability or other measure of membership in a particular class; or service job prioritization in which a service center receives service jobs and a scheduler ranks the jobs in order of importance based on a suitable prioritization criterion.

Each such application can be generally described as having a judge (e.g., the search engine, or the classifier, or the scheduler) that assigns a score (relevance ranking, or soft classification probability, or job priority) to items (web sites, or classes of a set of classes, or service jobs) based on a scoring criterion or algorithm (a search query processed by the search engine's relevance ranking basis, or an input object processed by a classifier algorithm, or a service job input to a scheduling algorithm). The obtained result is dependent upon the scoring criterion used by the judge, and in general different judges may use different scoring criteria and hence generate different rankings.

It is known to aggregate the rankings from different judges to produce a consensus aggregation. For example, meta-search engines have been deployed which input a received query to different Internet search engines, collect the rankings from the different Internet search engines, and apply an aggregation algorithm to generate a consensus ranking. Some Internet users find that meta-search engines provide better results than any individual search engine operating alone. Similarly, “meta-classifiers” combine classification results from different classifiers to generate a consensus classification.

In such consensus aggregation, an aggregation function is applied to generate consensus scores for the items based on the scores assigned by the constituent judges. The consensus ranking is dependent upon the choice of aggregation function. For example, a simple aggregation function is an average, in which the scores generated by the different judges for each item are averaged together to generate the aggregated score for that item. A disadvantage of this approach is that it can produce a high aggregated score for an item on which the judges have little consensus. For example, an item for which half the judges rank near the top of the ranking and half the judges rank near the bottom of the ranking will end up near the middle of the aggregated ranking. However, that is not reflective the consensus of the judges.

Some more complex aggregation functions weight the average score of an item by the number of judges that rank the item above a threshold (for example, the number of judges that rank the item in the “top ten”). These approaches are deemed to improve the aggregated rankings. However, the basis for this modification of the aggregation function is not readily apparent, and the approach can overemphasize a generally high ranking (e.g., in the “top ten”) compared with the actual scores assigned by the various judges. These aggregation functions also do not recognize or take into account relationships which may or may not exist between different judges, and are not tunable to accommodate different relationships between judges. For example, consensus between two judges that use wholly independent ranking criteria may be of substantial significance; whereas, consensus between two judges that use highly correlated ranking criteria may be of little significance.

BRIEF DESCRIPTION

In some illustrative embodiments disclosed as illustrative examples herein, a fusion method is disclosed for fusing M rankings generated by M judges where M is an integer greater than unity, the fusion method comprising: computing values of an aggregation function for items in a union of items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item; constructing an aggregation ranking of items in the union of items of the M rankings based on the computed values of the aggregation function; and outputting informational content embodying the aggregation ranking.

In some illustrative embodiments disclosed as illustrative examples herein, a storage medium is disclosed, which stores instructions executable to fuse M rankings generated by M judges where M is an integer greater than unity by (i) computing values of an aggregation function for items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item and (ii) constructing an aggregation ranking of items of the M rankings based on the computed values of the aggregation function.

In some illustrative embodiments disclosed as illustrative examples herein, a fusion apparatus is disclosed for fusing M rankings generated by M judges where M is an integer greater than unity, the fusion apparatus comprising a digital processing device configured to (i) compute values of an aggregation function for items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item and (ii) construct an aggregation ranking of items of the M rankings based on the computed values of the aggregation function.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 diagrammatically shows a fusion apparatus for fusing M rankings generated by M judges where M is an integer greater than unity.

FIG. 2 diagrammatically shows an application of the fusion apparatus of FIG. 1 to generating a consensus search result from search results generated by a query input to M different Internet search engines.

FIG. 3 diagrammatically shows an application of the fusion apparatus of FIG. 1 to generating a consensus probability vector from M probability vectors each having k elements generated for an input object by M different soft classifiers.

DETAILED DESCRIPTION

With reference to FIG. 1, a fusion apparatus is illustrated. The fusion apparatus is suitably embodied by a digital processing device 4, such as a computer (illustrated), or a personal data assistant (PDA), cellular telephone or other handheld device, or so forth. In some embodiments the digital processing device is Internet-based, for example embodied as an Internet server accessible via the Internet or an organizational server accessible via an organizational digital network. The fusion apparatus optionally also includes one or more output devices 6, such as a display (illustrated), a printer or other marking engine, a voice synthesizer, or so forth. The fusion apparatus optionally also includes one or more user input devices 8, such as a keyboard (illustrated), or a mouse, trackball, or other pointing-type input device, or so forth. The illustrated input and output devices 6, 8 suitably define a user interface. In the illustrated embodiment, the unitary computer 4 includes the input and output devices 6, 8 such that the processing device and the user interface device are integrated together. Alternatively, they may be separate—for example, the user interface device may be embodied as a laptop computer while the processing device performing the fusion operations is embodied as a remote Internet-based server system. The fusion system may also be embodied as a storage medium storing instructions executable to perform the disclosed fusion operations. For example, the storage medium may be a magnetic disk or tape, an optical disk or tape, an electrostatic memory, a random access memory (RAM), a read-only memory (ROM), a remote Internet server storage device, or so forth, and the stored instructions may be executable on one or more of the described digital processing devices such as a computer, PDA, cellular telephone, Internet server, or so forth.

The fusion apparatus communicates with M judges 10, where M is an integer greater than unity, to receive M raw rankings 12. For example, the judges 10 may be different Internet search engines and the raw rankings 12 may be search results generated by the search engines responsive to a query. In other applications, the judges 10 may be different soft classifiers each of which is configured to receive an input object and to output a probability vector of length or dimension k where the values of the k elements provide membership probabilities for the input object in each of k classes of a set of k classes. In general, each judge generates its ranking using its own criterion or criteria, which may be wholly different from the criterion or criteria used by other judges, or which may partially or wholly overlap the criterion or criteria used by one or more other judges. As a consequence, the M raw rankings 12 are generally different, although a given item may be ranked highly or at the same ranking value in two or more of the M raw rankings 12. For example, if the M judges are different Internet search engines and they receive a query “George Washington”, it may not be surprising that some, many, or all of the M Internet search engines may rank a site entitled “The life and times of George Washington” that is maintained by a respected historical organization near the top of their respective raw rankings 12. On the other hand, the rankings will differ to the extent that different criteria are used by different judges. To use the same example, if the site entitled “The life and times of George Washington” is a scholarly site that is infrequently visited except by a small and select group of scholars, then an Internet search engine that strongly emphasizes number of visits in its ranking criterion may rank this site rather low, whereas another search engines that place less emphasis on number of visits may rank the site substantially higher.

Depending upon the nature of the raw rankings 12, some preprocessing may optionally be performed in order to convert the raw rankings 12 into a format more amenable to aggregation. For example, a ranking truncator 14 may truncate each of the raw rankings 12 so as to retain only the top k results. Keeping with the search engines example, a given search engine may report dozens, hundreds, thousands, or more results, and in such a case the rankings truncator 14 may retain, for example, the top k=10 results in each raw ranking 12. On the other hand, if the judges are different classifiers each of which is trained to perform soft classification of objects respective to a set of C classes, then the raw rankings are suitably probability vectors each having precisely C elements whose values are indicative of membership probability in the corresponding C classes. In this case, truncation is optionally omitted (in which case k=C), or truncation may be performed to reduce the size of the rankings to a value k that is smaller than C. In the following, it is assumed that either: (i) the truncator 14 ensures that the truncated rankings each include k items; or (ii) the raw rankings 12 inherently each have k items. However, this assumption is only for convenience, and the described fusion techniques can also be applied to fuse rankings having different numbers of items.

The following notation is used herein. The set of all ranked items that are ranked in the top k by at least one of the M judges 10 is denoted as the set U, which is the union of the rankings (after truncation, if performed). The size of the set U is K, where K≧k holds. The condition K=k holds if every one of the rankings (after truncation if performed) contains the same k items (albeit possibly in different orders). The condition K>k holds if some items are ranked in one or more rankings but not in one or more other rankings (again, after truncation if performed). The rankings 12 are denoted by the notation L_(k)(j) where k denotes the ranking either inherently includes k items or has been truncated to include only the top k items, and j indexes the j^(th) judge.

As a further preprocessing consideration, different judges, such as different Internet search engines, may use different scales in ranking results. For example one search engine may rank a result on a scale of one-to-five stars, whereas another search engine may rank a result on a scale of 0-100%. Accordingly, in such applications a ranking normalizer 16 operates to normalize the rankings to a common ranking scale. A suitable normalization is the score normalization, which is suitably implemented as follows. Defining V_(j)(i) as the unnormalized ranking value assigned by the j^(th) judge to the i^(th) item of the set of set of ranked items U, the normalized ranking value S_(j)(i) is:

$\begin{matrix} {{{S_{j}(i)} = \frac{{V_{j}(i)} - {\min_{i^{\prime}}\left\{ {V_{j}\left( i^{\prime} \right)} \right\}}}{{\max_{i^{\prime}}\left\{ {V_{j}\left( i^{\prime} \right)} \right\}} - {\min_{i^{\prime}}\left\{ {V_{j}\left( i^{\prime} \right)} \right\}}}}{{{\forall j} = 1},\ldots \mspace{11mu},{M;{\forall{i \in {\mathcal{L}_{k}(j)}}}},}} & (1) \end{matrix}$

for items that are included in the ranking L_(k)(j), and:

S _(j)(i)=0 ∀j=1, . . . ,M; ∀iε(U−L _(k)(j)   (2)

for items that are not included in the ranking L_(k)(j). Equation (1 ) ensures that the ranked values S_(j)(i) lie in the interval [0,1] for items that are included in the ranking L_(k)(j), and Equation (2) ensures that S_(j)(i)=0 for items that are not included in the ranking L_(k)(j) (but which are ranked in at least one other ranking L_(k)(j′) where j′≠j, and hence are members of the set of ranked items U).

The rankings after appropriate preprocessing such as truncation and/or normalization are shown in FIG. 1 as the set of M rankings 20. The union of these M rankings 20 is the set of ranked items U 22, which includes K items. The rankings 20 have values suitably denoted as S_(j)(i) for the i^(th) item (in the list of ranked items U 22), where any item that does not appear in the j^(th) ranking (and hence is deemed “not relevant” by judge j) is assigned a value S_(j)(i)=0 by operation of Equation (2).

The M rankings 20 and the set of ranked items U 22 serve as input to an aggregation function 30 that is used by an items ranker 32 to generate a consensus or aggregation ranking 34 of the set of ranked items U 22. In a suitable embodiment, the aggregation function 30 operates as follows. The parameter S_(j)(i) as used herein denotes an event in which the j^(th) judge considers the i^(th) item to be relevant, that is, included in the ranking L_(k)(j). The normalized score S_(j)(i) has a value in the range [0,1], and can be treated as a truth value of event S_(j)(i). The parameter E_(m) ^(M)(i) herein denotes an event in which at least m judges consider the i^(th) item to be relevant, that is, in which the i^(th) item is included in at least m of the rankings L_(k)(i), j=1, . . . ,M. The event E_(m) ^(M)(i) is given by:

$\begin{matrix} {{E_{m}^{M} = {\bigcup\limits_{1 \leq j_{1} < j_{2} < \; \ldots \mspace{11mu} < j_{m} \leq M}\left( {\bigcap\left( {{S_{j_{1}}(i)},{S_{j_{2}}(i)},,\ldots \mspace{14mu},{S_{j_{m}}(i)}} \right)} \right)}}{{{\forall m} = 1},\ldots \mspace{11mu},{M;{{\forall i} = 1}},\ldots \mspace{11mu},{K.}}} & (3) \end{matrix}$

The truth value of event E_(m) ^(M)(i) is denoted herein as E_(m) ^(M)(i). This is a combinatorial problem which generalizes Poincaré's formula also known as the inclusion-exclusion principle. It can be shown that the combinatorial computation of E_(m) ^(M)(i) can be expressed as a linear combination of the following quantities, ∀i=1, . . . ,K:

$\begin{matrix} {{S_{l}^{M}(i)} = \left\{ \begin{matrix} {\sum\limits_{j = 1}^{M}{S_{j}(i)}} & {{{for}\mspace{14mu} l} = 1} \\ {\sum\limits_{\underset{< j_{l} \geq M}{1 \leq j_{1} < j_{2} < \mspace{11mu} \ldots}}\left( {T\left( {{S_{j_{1}}(i)},{S_{j_{2}}(i)},\ldots \mspace{14mu},{S_{j_{l}}(i)}} \right)} \right)} & {{{{for}\mspace{14mu} l} = 2},\ldots \mspace{14mu},M} \end{matrix} \right.} & (4) \end{matrix}$

where T is in some embodiments a triangular norm, and more generally is a conjunction operation.

Triangular norms are logical operations which, given two events A and B with truth values μ(A) and μ(B) both in the range [0, 1], give a truth value, also in the range [0, 1], of the conjunction A∩B. Triangular norms are also sometimes called t-norms, and are special functions from [0, 1]×[0, 1] to [0, 1] which exhibit certain characteristics such as being commutative, associative, and monotonic.

The linear combination of Equation (4) is related to Jordan's combinatorial formulas, and can be written as:

$\begin{matrix} {{{E_{m}^{M}(i)} = {\sum\limits_{l = m}^{M}{\left( {- 1} \right)^{l - m}\begin{pmatrix} {l - 1} \\ {m - 1} \end{pmatrix}{S_{l}^{M}(i)}}}}{{{\forall m} = 1},\ldots \mspace{14mu},{M;}}{{{\forall i} = 1},\ldots \mspace{14mu},K}} & (5) \end{matrix}$

where expressions of the form:

$\begin{matrix} {\begin{pmatrix} n \\ p \end{pmatrix} = \frac{n!}{{p!}{\left( {n - p} \right)!}}} & (6) \end{matrix}$

are binomial coefficients. Equation (5) can be further rewritten as:

$\begin{matrix} {{{E_{m}^{M}(i)} = {\sum\limits_{l = m}^{M}{\begin{pmatrix} {- m} \\ {l - m} \end{pmatrix}{S_{l}^{M}(i)}}}}{{{\forall m} = 1},\ldots \mspace{14mu},{M;{{\forall i} = 1}},\ldots \mspace{14mu},K}} & (7) \end{matrix}$

where expressions of the form:

$\begin{matrix} {\begin{pmatrix} {- n} \\ p \end{pmatrix} = {{\left( {- 1} \right)^{p}\begin{pmatrix} {n + p - 1} \\ p \end{pmatrix}} = \frac{\left( {- n} \right)\left( {{- n} - 1} \right)\mspace{20mu} \ldots \mspace{14mu} \left( {{- n} - p + 1} \right)}{p!}}} & (8) \end{matrix}$

are negative binomial coefficients.

The truth values E_(m) ^(M)(i); m=1, . . . ,M can be viewed as consensus measures indicative of the consensus amongst the M judges that the i^(th) item should be ranked as relevant. A consensus aggregation function A(i) is constructed from these consensus measures E_(m) ^(M)(i);m=1, . . . ,M as follows:

$\begin{matrix} {{{A(i)} = {{\sum\limits_{m = 1}^{M}{{{mE}_{m}^{M}(i)}\mspace{20mu} {\forall i}}} = 1}},\ldots \mspace{11mu},{K.}} & (9) \end{matrix}$

It is desired that the i^(th) item should have a higher aggregation score as the number of judges that find the i^(th) item to be relevant increases. In accordance with this desire, the aggregation function of Equation (9) weighs by m the consensus measure E_(m) ^(M)(i) that item i is relevant for at least m judges, and takes a sum of these weighted values over m from 1 to M. This approach gives increasing weights with increasing value of m, that is, with increasing number of judges who consider item i to be relevant.

The computational cost of Equation (9) is relatively high. Computation of S_(l) ^(M)(i) in Equation (4) entails enumeration of

$\left. \left( \begin{matrix} M \\ l \end{matrix}\quad \right. \right)$

combinations which grows exponentially. However, it can be shown that Equation (9) can be rewritten as:

$\begin{matrix} {{A(i)} = {{\sum\limits_{j = 1}^{M}{S_{j}(i)}} + {\sum\limits_{1 \leq j < j^{\prime} \leq M}{T\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)}}}} & (10) \end{matrix}$

which has substantially lower computational cost of order O(M²). In Equation (10), T(. . .) denotes a pairwise conjunction operator 40. In some embodiments described herein, T(. . .) denotes a triangular norm conjunction operator.

In some suitable embodiments, the pairwise conjunction operator T(. . .) 40 is a non-parametric triangular norm, such as for example: a minimum triangular norm T^(M) given by:

T ^(M)(S _(j)(i),S _(j′)(i)=min(S _(j)(i), S _(j′)(i))   (11);

a product triangular norm T^(P) given by:

T ^(P)(S _(j)(i),S _(j′)(i))=(S _(j)(i)·S _(j′)(i)) where · denotes a product   (12);

a Lukasiewicz triangular norm T^(L) given by:

T ^(L)(S _(j)(i),S _(j′)(i)=max(S _(j)(i)+S _(j′)(i)−1, 0)   (13);

or a drastic product triangular norm T^(D) given by:

$\begin{matrix} {{T^{D}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} = \left\{ {\begin{matrix} 0 & {{{if}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} \in \left\lbrack {0,{1\left\lbrack {}^{2} \right.}} \right.} \\ {\min \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} & {otherwise} \end{matrix}.} \right.} & (14) \end{matrix}$

Employing a product triangular norm T^(P) for the pairwise conjunction of ranking values of different judges has the effect of assuming that the pair of judges are independent, that is, have a low or zero strength of relationship. On the other hand, employing a minimum triangular norm T^(M) for the pairwise conjunction of ranking values of different judges has the effect of assuming that the pair of judges are highly correlated, that is, have a high strength of relationship. In probability theory, T^(M) and T^(L) are respectively upper and lower bounds for copulas. They are called Fréchet-Hoeffding boundaries. In triangular norms theory, T^(M) is the upper bound of triangular norms but T^(D) is the one which is the lower bound. That is:

T ^(D)(μ(A),μ(B)≦T(μ(A),μ(B)≦T ^(M)(μ(A),μ(B)) ∀T(. . .) is a triangular norm; ∀(μ(A),μ(B))ε[0,1[²  (15);

As T^(M) is the upper bound of triangular norms, this means that it gives the highest conjunction measure between two judges. It can be shown that by using T^(M) as the relationship type between two judges j, j′, it follows that either S_(j)(i)⊂S_(j′)(i) or S_(j′)(i)⊂S_(j)(i). On the other hand, as T^(D) is the lower bound of triangular norms, it typically yields a low or zero value when used as the pairwise conjunction.

Using non-parametric t-norms for the pairwise conjunctions can reduce the computation time. For example, the following factorization property holds:

$\begin{matrix} {{\sum\limits_{1 \leq j < j^{\prime} \leq M}{T^{P}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)}} = {\frac{1}{2}\left( {\left( {\sum\limits_{j = 1}^{M}{S_{j}(i)}} \right)^{2} - {\sum\limits_{j = 1}^{M}\left( {S_{j}(i)} \right)^{2}}} \right)}} & (16) \end{matrix}$

which leads to O(M) computation cost when using the product triangular norm T^(P). The computational cost when using the minimum triangular norm T^(M) can be reduced to O(M log(M)) by using the following:

$\begin{matrix} {{{\sum\limits_{1 \leq j < j^{\prime} \leq M}{T^{M}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)}} = {\sum\limits_{j = 1}^{M}{\left( {M - j} \right){S_{j}^{*}(i)}}}}{where}{{S_{1}^{*}(i)} \leq {S_{2}^{*}(i)} \leq \ldots \leq {{S_{M}^{*}(i)}.}}} & (17) \end{matrix}$

In some suitable embodiments, the pairwise conjunction operator T(. . .) 40 of Equation (10) is a parametric triangular norm, such as for example: a Schweizer-Sklar parametric triangular norm defined for the parameter range ∀λε[−∞,∞] and given by

$\begin{matrix} {{T_{\lambda}^{SS}\left( {{\mu (A)},{\mu (B)}} \right)} = \left\{ {\begin{matrix} {T^{M}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = {- \infty}} \\ {T^{P}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = 0} \\ {T^{D}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = \infty} \\ \left( {\max \left( {{{\mu (A)}^{\lambda} + {\mu (B)}^{\lambda} - 1},0} \right)} \right)^{\frac{1}{\lambda}} & {otherwise} \end{matrix};} \right.} & (18) \end{matrix}$

a Frank triangular norm defined for the parameter range ∀λε[0,∞] and given by

$\begin{matrix} {{T_{\lambda}^{F}\left( {{\mu (A)},{\mu (B)}} \right)} = \left\{ {\begin{matrix} {T^{M}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = 0} \\ {T^{P}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = 1} \\ {T^{L}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = \infty} \\ {\log_{\lambda}\left( {1 + \frac{\left( {\lambda^{\mu {(A)}} - 1} \right)\left( {\lambda^{\mu {(B)}} - 1} \right)}{\lambda - 1}} \right)} & {otherwise} \end{matrix};} \right.} & (19) \end{matrix}$

or an Aczél-Alsina triangular norm defined for the parameter range ∀λε[0,∞] and given by

$\begin{matrix} {{T_{\lambda}^{AA}\left( {{\mu (A)},{\mu (B)}} \right)} = \left\{ {\begin{matrix} {T^{D}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = 0} \\ {T^{M}\left( {{\mu (A)},{\mu (B)}} \right)} & {{{if}\mspace{14mu} \lambda} = \infty} \\ {\exp \left( {- \left( {\left( {{- \log}\; {\mu (A)}} \right)^{\lambda}\left( {{- \log}\; {\mu (B)}} \right)^{\lambda}} \right)^{1/\lambda}} \right)} & {otherwise} \end{matrix}.} \right.} & (20) \end{matrix}$

These parametric triangular norms reduce to one of the non-parametric triangular norms of Equations (11)-(14) for certain values of the parameter λ. For example, T_(λ=0) ^(AA)(μ(A),μ(B))=T^(D)(μ(A),μ(B)) for λ=0.

The aggregation function A(i) as written in Equation (10) expresses A(i) in terms of pairwise conjunctions between judges. This optionally enables specification of different types of the relationships between different pairs of judges by using different fuzzy conjunction operations. For example, the choice of fuzzy conjunction operation can be made to reflect a strength of relationship between the pair of judges. When using a parametric triangular norm such as one of the illustrative parametric triangular norms of Equations (18)-(20), this can be achieved by having a strength of relationship estimator 42 estimate strengths of relationship for different pairs of ranking values obtained from the M rankings, and adjusting a parameter of the parametric triangular norm (λ, in the illustrative examples of Equations (18)-(20)) for different pairwise conjunctions of ranking values based on the estimated strengths of relationship of the different pairs of ranking values.

For example, when using parametric copulas such as T_(λ) ^(AA), the adjustment entails providing different pairwise parameters λ_(jj′);∀j<j′ to reflect the different strengths of relationships between different pairs of judges. As the parameters λ_(jj′) give the strength of the relationship between judges j and j′, these values are optionally computed using Kendall correlation measures. In that configuration, each S_(j);j=1, . . . ,M is interpreted as a random variable and for a pair of judges j, j′ their different profiles are compared by computing their Kendall correlation measure given by:

$\begin{matrix} {{\tau \left( {S_{j},S_{j^{\prime}}} \right)} = \frac{{{Conc}\left( {S_{j},S_{j^{\prime}}} \right)} - {{Disc}\left( {S_{j},S_{j^{\prime}}} \right)}}{\sqrt{\left( {\frac{N\left( {N - 1} \right)}{2} - {{Tie}\left( S_{j} \right)}} \right)\left( {\frac{N\left( {N - 1} \right)}{2} - {{Tie}\left( S_{j^{\prime}} \right)}} \right)}}} & (21) \\ {{{{where}:{{Conc}\left( {S_{j},S_{j^{\prime}}} \right)}} = {\# \left\{ {\left( {i,i^{\prime}} \right),{i < {i^{\prime}:{\left\{ {{S_{j}(i)} < {{S_{j}\left( i^{\prime} \right)}\bigcap{S_{j^{\prime}}(i)}} < {S_{j^{\prime}}\left( i^{\prime} \right)}} \right\}\bigcup\left\{ {{S_{j}(i)} > {{S_{j}\left( i^{\prime} \right)}\bigcap{S_{j^{\prime}}(i)}} > {S_{j^{\prime}}\left( i^{\prime} \right)}} \right\}}}}} \right\}}},} & (22) \\ {{{{Disc}\left( {S_{j},S_{j^{\prime}}} \right)} = {\# \left\{ {\left( {i,i^{\prime}} \right),{i < {i^{\prime}:{\left\{ {{S_{j}(i)} < {{S_{j}\left( i^{\prime} \right)}\bigcap{S_{j^{\prime}}(i)}} > {S_{j^{\prime}}\left( i^{\prime} \right)}} \right\}\bigcup\left\{ {{S_{j}(i)} > {{S_{j}\left( i^{\prime} \right)}\bigcap{S_{j^{\prime}}(i)}} < {S_{j^{\prime}}\left( i^{\prime} \right)}} \right\}}}}} \right\}}},} & (23) \\ {{and}{{{Tie}\left( {S_{j},S_{j^{\prime}}} \right)} = {\# \left\{ {\left( {i,i^{\prime}} \right),{i < {i^{\prime}:{\left\{ {{S_{j}(i)} = {S_{j}\left( i^{\prime} \right)}} \right\}.}}}} \right.}}} & (24) \end{matrix}$

In probability theory, a copula is a function that takes as inputs the marginal distribution of two random variables and gives as output their joint distribution. A copula characterizes the dependence type between two random variables. In that context, for some parametric copulas a link function exists that allows computation of the parameter λ_(jj′) from the Kendall correlation measure τ(S_(j),S_(j′)). As an illustrative example, the strength of relationship estimator 42 suitably computes the following link function for the Aczél-Alsina triangular norms T_(λ) ^(AA) of Equation (20):

$\begin{matrix} {\lambda_{{jj}^{\prime}} = {\frac{1}{1 - {\tau \left( {S_{j},S_{j^{\prime}}} \right)}}.}} & (25) \end{matrix}$

Estimating the parameter λ_(jj′) by computing pairwise Kendall correlation measures is not mathematically rigorous for triangular norms, but is suitable for computing automatic adapted values for each pair λ_(jj′);∀j<j′.

It is also contemplated for the strength of relationship estimator 42 to use different correlation measures for estimating the relationships between the profiles of different pairs of judges. It is advantageous to have, for each judge j, a distribution of correlation measures which identifies the judge or judges j′ which are the most similar with judge j, in order to give more emphasis to their joint measure. Toward this end, either a Bravais-Pearson correlation measure or a Spearman's rank correlation measure can be associated with an adapted link function for computing the parameters λ_(jj′).

In the foregoing, various triangular norm functions have been described which are suitable for use as the pairwise conjunction operator T(. . .) 40 in the consensus aggregation function A(i) as set forth in Equation (10). However, substantially any kind of binary conjunction operation that measures the conjunction of two fuzzy sets can also be used as the pairwise conjunction operator T(. . .) 40 in Equation (10). As an example, the prioritized conjunction of Dubois and Prade:

F ^(DP)(S _(j)(i),S _(j′)(i)=min(S _(j)(i),max(S _(j′)(i), 1−h(S _(j) ,S _(j′))))  (26)

can be used as the pairwise conjunction operator T(. . .) 40, where

${h\left( {S_{j},S_{j^{\prime}}} \right)} = {\sup\limits_{i}\left( {\min \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} \right)}$

is the degree of consistency between two judges.

The fusion apparatus of FIG. 1 for fusing M rankings generated by M judges can be used in various ways. In the following, two illustrative applications are illustrated and described.

With reference to FIG. 2, an illustrative meta-search engine application is illustrated. In this application, a search query 100 is input to M Internet search engines 110 which serve as the M judges 10 of FIG. 1. The M Internet search engines 110 generate M respective search results 112, which serve as the M raw rankings 12 of FIG. 1. The fusion apparatus 14, 16, 20, 22, 30, 32 operates on the M search results 112 as described with reference to FIG. 1 to generate a consensus search result 134 that corresponds to the aggregation ranking 34 of FIG. 1.

In one suitable configuration, the M Internet search engines 110 are publicly accessible search engines such as Google™, Altavista™, Lycos™, or so forth, and the fusion apparatus forwards the query 100 to each search engine at its universal resource locator (URL) address (for example, http://www.google.com/ for Google™) using a query string of the form “?<query>” appended to the URL, or using another suitable syntax, to input the query 100 to the search engine. In some embodiments, the fusion apparatus 14, 16, 20, 22, 30, 32 is embodied as an Internet server that is accessed by a client computer using a suitable URL and provides a user interface configured to receive the query 100 from the client computer via the Internet. The Internet server embodying the fusion apparatus 14, 16, 20, 22, 30, 32 then computes the consensus search result 134 and returns the consensus search result 134 to the client computer via the Internet in a suitable format such as HTML.

In other contemplated embodiments, the fusion apparatus 14, 16, 20, 22, 30, 32 may be a plug-in or other component of a web browser executing on a laptop computer, desktop computer, or other computer directly accessed by an end user. The plug-in or other component of the web browser is configured to receive the query 100 from the user, input the query 100 to each of the M Internet search engines 110 over the Internet via the respective URL addresses with a query string of the form “?<query>” appended to each URL, or using another suitable Internet syntax, and receives the M search results 112 from the Internet search engines 110. Thereafter, the fusion apparatus 14, 16, 20, 22, 30, 32 embodied as a plug-in or other component of the web browser computes and displays the consensus search result 134.

With reference to FIG. 3, in a soft classification application example an input object 200 is input to M soft classifiers 210 which in this application serve as the M judges 10 of the fusion apparatus of FIG. 1. The M soft classifiers 210 generate M respective probability vectors 212 corresponding to the M raw rankings 12 of FIG. 1. In this illustrative embodiment, the each of the M soft classifiers 210 are configured to classify the input object 200 respective to classes of a set of k classes, and each of the M soft classifiers 210 is further configured to generate a normalized probability vector in which the k elements of the probability vector are each in the range [0,1] and the sum of the k elements of the probability vector is unity. Accordingly, in this embodiment the truncator 14 and the ranking normalizer 16 are suitably omitted, since each of the M probability vectors 212 corresponding to the M raw rankings 12 of FIG. 1 are inherently of length k and are inherently normalized. The thusly simplified fusion apparatus 22, 30, 32 outputs a consensus probability vector 234 corresponding to the aggregation ranking 34 of FIG. 1.

In one suitable embodiment, the M soft classifiers 210 and the fusion apparatus 22, 30, 32 are all embodied by a single computer that runs software implementing the M soft classifiers 210 and the fusion apparatus 22, 30, 32 and also implementing a user interface that enables the user to identify or input the input object 200 and to perceive the consensus probability vector 234 in a human-perceptible format such as a list of displayed probabilities. In some embodiments, the software further implements hard classification by performing thresholding on the consensus probability vector 234 to assign the input object 200 to a single “most probable” class of the set of k classes. Moreover, it is contemplated for the M soft classifiers 210 and the fusion apparatus 22, 30, 32 to be an automated composite component of larger overall system such that the fusion apparatus 22, 30, 32 does not directly interact with a user. For example, the M soft classifiers 210 and the fusion apparatus 22, 30, 32 may be components of a document storage and retrieval system, the input object may be an input document to be stored, and the consensus probability vector 234 is not displayed but rather is used by the document storage and retrieval system to (i) index or otherwise organize the input document respective to other documents stored in the document storage and retrieval system, and to (ii) retrieve the input document at a later date in response to receipt of a user query criterion that is satisfied by the consensus probability vector 234.

It is to be appreciated that FIGS. 2 and 3 merely illustrate some applications as examples. The fusion apparatus as described with reference to FIG. 1 will find use in any application in which a plurality of judges rank items by order of relevance, such that the different judges produce different rankings from which it is advantageous to generate an aggregate or consensus ranking.

It will be appreciated that various of the above-disclosed and other features and functions, or alternatives thereof, may be desirably combined into many other different systems or applications. Also that various presently unforeseen or unanticipated alternatives, modifications, variations or improvements therein may be subsequently made by those skilled in the art which are also intended to be encompassed by the following claims. 

1. A fusion method for fusing M rankings generated by M judges where M is an integer greater than unity, the fusion method comprising: computing values of an aggregation function for items in a union of items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item; constructing an aggregation ranking of items in the union of items of the M rankings based on the computed values of the aggregation function; and outputting informational content embodying the aggregation ranking.
 2. The fusion method as set forth in claim 1, further comprising: prior to the computing, truncating the M rankings to include only k highest ranked items of each truncated ranking where k is an integer greater than unity; and normalizing the M truncated rankings.
 3. The fusion method as set forth in claim 1, wherein the pairwise conjunctions are pairwise triangular norms.
 4. The fusion method as set forth in claim 1, wherein the aggregation function has an output quantitatively equal to A(i) where A(i) is given by: ${A(i)} = {{\sum\limits_{j = 1}^{M}{S_{j}(i)}} + {\sum\limits_{1 \leq j < j^{\prime} \leq M}{T\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)}}}$ where the index i indexes the input item in the union of items of the M rankings, S_(j)(i) denotes the ranking value for input item i generated by a judge indexed by j of the M judges, and T(. . .) denotes a triangular norm.
 5. The fusion method as set forth in claim 1, wherein the sum of pairwise conjunctions of ranking values of different judges for an input item included in the aggregation function has an output quantitatively equal to an output of the sum: $\sum\limits_{1 \leq j < j^{\prime} \leq M}{T\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)}$ where the index i indexes the input item in the union of items of the M rankings, S_(j)(i) denotes the ranking value for input item i generated by a judge indexed by j of the M judges, and T(. . .) denotes a triangular norm.
 6. The fusion method as set forth in claim 5, wherein the triangular norm T(S_(j)(i),S_(j′)(i)) is a non-parametric triangular norm selected from a group consisting of: min(S_(j)(i),S_(j′(i)),) (S_(j)(i)·S_(j′)(i)) where · denotes a product, a Lukasiewicz norm T^(L)(S_(j)(i),S_(j′)(i))=max(S_(j)(i)+S_(j′)(i)−1, 0) and a drastic product ${T^{D}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} = \left\{ {\begin{matrix} 0 & {{{if}\mspace{14mu} \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} \in \left\lbrack {0,{1\left\lbrack {}^{2} \right.}} \right.} \\ {\min \; \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} & {otherwise} \end{matrix}.} \right.$
 7. The fusion method as set forth in claim 5, wherein the triangular norm T(S_(j)(i),S_(j′)(i)) is a parametric triangular norm.
 8. The fusion method as set forth in claim 7, further comprising: estimating strengths of relationship for different pairs of ranking values obtained from the M rankings; and adjusting a parameter of the parametric triangular norm for different pairwise conjunctions of ranking values based on the estimated strengths of relationship of the different pairs of ranking values.
 9. The fusion method as set forth in claim 1, further comprising: estimating strengths of relationship for different pairs of ranking values for items in the union of items of the M rankings; and adjusting the pairwise conjunctions of ranking values of the sum of pairwise conjunctions of ranking values included in the aggregation function based on the estimated strengths of relationship for the different pairs of ranking values.
 10. The fusion method as set forth in claim 1, wherein the M judges are M different Internet search engines and the M rankings are M sets of search engine results generated for a query input to the M different Internet search engines, and the outputting of informational content embodying the aggregation ranking comprises: displaying a consensus Internet search result defined by the aggregation ranking and corresponding to the query.
 11. The fusion method as set forth in claim 1, wherein the M judges are M different soft classifiers for classifying an object respective to a set of k classes where k is an integer greater than unity, and the M rankings are M probability vectors each having k elements generated for an input object by the M different soft classifiers, and the outputting of informational content embodying the aggregation ranking comprises: assigning a classification to the input object based on a consensus probability vector defined by the aggregation ranking; and at least one of displaying and storing the classification.
 12. A storage medium storing instructions executable to fuse M rankings generated by M judges where M is an integer greater than unity by (i) computing values of an aggregation function for items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item and (ii) constructing an aggregation ranking of items of the M rankings based on the computed values of the aggregation function.
 13. The storage medium as set forth in claim 12, wherein the stored instructions are further executable to truncate the M rankings to include only k highest ranked items of each truncated ranking where k is an integer greater than unity and to normalize the M truncated rankings, the computing (i) and constructing (ii) operating on the M truncated rankings.
 14. The storage medium as set forth in claim 12, wherein the pairwise conjunctions are pairwise triangular norms.
 15. The storage medium as set forth in claim 14, wherein the pairwise triangular norms are selected from a group consisting of: ${T^{M} = {\min \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)}},{T^{P} = {\left( {{S_{j}(i)} \cdot {S_{j^{\prime}}(i)}} \right)\mspace{14mu} {{where} \cdot {denotes}}\mspace{14mu} a\mspace{14mu} {product}}},{{T^{L}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} = {\max \left( {{{S_{j}(i)} + {S_{j^{\prime}}(i)} - 1},0} \right)}},{{T^{D}\left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} = \left\{ {\begin{matrix} 0 & {{{if}\mspace{14mu} \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} \in \left\lbrack {0,{1\left\lbrack {}^{2}\mspace{14mu} \right.}} \right.} \\ {\min \left( {{S_{j}(i)},{S_{j^{\prime}}(i)}} \right)} & {otherwise} \end{matrix},{and}} \right.}$ a parametric triangular norm incorporating at least one of T^(M),T^(P),T^(L),T^(D), where index i indexes the input item and S_(j)(i) denotes the ranking value for input item i generated by a judge indexed by j of the M judges.
 16. The storage medium as set forth in claim 12, wherein the computing of values of the aggregation function for items of the M rankings further includes (i)(a) estimating strengths of relationship for different pairs of ranking values for items of the M rankings and (i)(b) adjusting the pairwise conjunctions of ranking values of the sum of pairwise conjunctions based on the estimated strengths of relationship for the different pairs of ranking values.
 17. The storage medium as set forth in claim 12, wherein the stored instructions are further executable to acquire the M rankings as M sets of search engine results generated for a query input to M different Internet search engines.
 18. The storage medium as set forth in claim 12, wherein the stored instructions are further executable to acquire the M rankings as M probability vectors each having k elements corresponding to classes of a set of k classes generated for an input object by M different soft classifiers.
 19. A fusion apparatus for fusing M rankings generated by M judges where M is an integer greater than unity, the fusion apparatus comprising: a digital processing device configured to (i) compute values of an aggregation function for items of the M rankings, the aggregation function including a sum of pairwise conjunctions of ranking values of different judges for an input item and (ii) construct an aggregation ranking of items of the M rankings based on the computed values of the aggregation function.
 20. The fusion apparatus as set forth in claim 19, further comprising: a user interface device integrated with or separate from the digital processing device, the user interface device configured to display informational content embodying the aggregation ranking. 